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Abstract. The hierarchical clustering observed in cold dark 
matter simulations results in highly clumped galactic halos. 
If the dark matter in our halo is made of weakly interact- 
ing massive particles (WIMPs), their annihilation products 
should be detectable in the higher density and nearby clumps. 
We consider WIMPs to be neutralinos and calculate the synch- 
rotron flux from their annihilation products in the presence of 
the Galactic magnetic field. We derive a self-consistent emis- 
sion spectrum including pair annihilation, synchrotron self- 
absorption, and synchrotron self-Compton reactions. The re- 
sulting radiation spans microwave frequencies that can be ob- 
served over the anisotropies in the cosmic microwave back- 
ground. These synchrotron sources should be identifiable as 
WIMP clumps, either by their spatial structure or by their 
distinctive radio spectrum. 



1 Introduction 

The density of dark matter in the Universe is observed via its 
gravitational effects on galaxies and clusters of galaxies to 
constitute about 30% of the critical density of the Universe. 
The nature of this dominant matter component is still un- 
known. Primordial nucleosynthesis and acoustic peaks in the 
cosmic microwave background (CMB) constrain the density 
of baryonic matter to be less than about 5% of the critical 
density, thus most of the dark matter is non-baryonic. The 
leading candidate for the dark matter is the lightest super- 
symmetric particle in supersymmetric extensions of the stan- 
dard model that can be stable by conservation of R-parity. 
In most scenarios this weakly interacting massive particle 
(WIMP) is the neutralino, which is a linear combination 
of the photino, the zino, and the higgsinos (for a review, see 
Jungman, Kamionkowski and Griest (1996)). 

Neutralinos may be detected directly as they traverse the 
Earth or indirectly by the observation of their annihilation 
products. Direct neutralino searches are now underway in 
a number of low temperature experiments with no consen- 
sus detection as of yet. Indirect searches have been pro- 
posed both for gamma rays and synchrotron emission from 
the annihilation of WIMPs in the Galactic center (Berezin- 
sky, Gurevich, and Zybin, 1992; Berezinsky, Bottino, and 
Mignola, 1994). The rate of neutralino annihilation is pro- 
portional to the neutralino density squared (~ n^.), therefore 
the strongest flux is expected to come from the highest den- 
sity regions such as the Galactic center. Depending on the 
history of formation of the central black hole in our Galaxy, 



the WIMP density may be strongly enhanced in the neigh- 
borhood of the black hole. The density enhancement cou- 
pled with assumptions about the local strength of the mag- 
netic field can, in principle, constrain viable models for neut- 
ralino masses and cross sections (Gondolo and Silk, 1999; 
Gondolo, 2000). 

Although promising, the proposal to focus indirect sear- 
ches on the Galactic nucleus is not free of uncertainties. The 
presence of central cusps in large halos has been questioned 
by observations Salucci (2000). In addition, the large frac- 
tion of baryonic matter as well as the history of the cen- 
tral black hole are likely to alter the dynamics of the dark 
matter in the Galactic center. Major mergers in the central 
black hole's past will tend to decrease the dark matter den- 
sity. A clear picture of the magnetic structure around a cen- 
tral black hole is also lacking. In addition, when the neutral- 
ino density becomes very high, the synchrotron signature is 
strongly modified by synchrotron self-Compton scatterings 
and pair annihilation, as we discuss below and in (Aloisio et 
al., 2002a). A cleaner dark matter annihilation signal may 
be detected against a better understood background using the 
clumpy nature of our halo. 

Superimposed on the smooth component of CDM halos, 
high-resolution simulations find a large degree of substruc- 
ture formed by the constant merging of smaller halos to form 
a present dark matter halo (see, for example, Ghigna et al. 
(1998)). The large number of clumps generated through the 
hierarchical clustering of CDM comprises about 10 — 20% 
of the total mass of a given halo. High density and nearby 
clumps in our own halo enhance appreciably the emission 
of gamma rays and neutrinos from neutralino annihilation 
(Bergstromet al., 1999; Calcaneo-Roldan and Moore, 2001). 
As these clumps reach closer to the Galactic plane, the Galac- 
tic magnetic field strength increases drastically and the synch- 
rotron signal from the charged products of neutralino annih- 
ilation intensifies considerably. 

Here, we calculate in detail the synchrotron radiation of 
electrons and positrons generated as decay products of WIMP 
annihilation in the Galactic magnetic field. We show that 
the synchrotron emission can provide a crucial test of WIMP 
models, since the predicted fluxes are in the microwave re- 
gion and exceed the signal of CMB anisotropies at some fre- 
quencies. The detection of this excess radiation from small 
angular size regions in the sky may provide the first signal of 
WIMP dark matter. 

In section 2, we explain the calculation procedure. In sec- 
tions 3 and 4 we present results for two different clump den- 



sity profiles. Section 5 is the conclusion, with prospects for 
detection. 



2 Synchrotron from Dark Matter Clumps 

2. 1 Dark Matter Clump Structure 

The density profile of dark matter halos has yet to be pre- 
cisely determined. The possible profiles for the smooth dark 
matter component on the scales of galaxies is reasonably 
well constrained by observations of galactic dynamics to- 
gether with CDM simulations. In contrast, the structure of 
the smaller dark matter halos (i.e., clumps) is only constrained 
by theoretical arguments and numerical simulations. Our 
Galactic halo can be described by both a smooth large scale 
component plus a distribution of clumps which are the smaller 
halos that have fallen into the Galactic potential well as it 
formed from lower mass objects. 

In order to derive the synchrotron emission from neutral- 
ino annihilation in the dark matter clumps of our Galaxy, 
we need to model both the smooth component of the halo 
as well as the individual density profile of the clumps. We 
consider the halo of the Galaxy to be well described by the 
Navarro, Frenk, and White (NFW) density profile (Navarro, 
Frenk and White, 1996, 1997) which we describe below. We 
model the individual clumps of dark matter in our halo by 
two choices of CDM halo profiles that span the range of rea- 
sonable slopes: the NFW profile and the singular isothermal 
sphere (SIS). These two models bracket a range of possible 
clump density profiles. 

The NFW profile arises from CDM simulations and is con- 
sistent with observations of a number of galaxies. The dark 
matter density as a function of radius from the center of the 
particular halo is given by: 



p(r) = m x n — 



1 + ^ 



(1) 



where r c is a core radius and m x is the neutralino mass. 
This profile gives a broad core region with a gently slop- 
ing profile inside (p ~ r _1 ) and a steeper slope outside r c 
(p ~ r~ 3 ). In general, n and r c are fixed by the total mass 
of the halo inside some radius i?haio 3> r c and the density at 
some particular point (or equivalently the velocity dispersion 
at that point). To set these parameters for the smooth com- 
ponent of our Galactic halo, we chose r c = i?haio/10 where 
-Rhaio is the radius within which lies most of our halo's mass, 
M halo (r < -Rhaio) = 2 x 10 12 M Q . In addition, we set the 
density in the solar neighborhood from the local velocity dis- 
persion to be p Q = p(8.5kpc) = 6.5 x 10~ 25 g/cm 3 which 
fixes all parameters. 

For each individual clump described by the NFW profile, 
we fix r c and n by setting r c = R c /10 where R c is the 
radius at which the density of the clump equals that of the 
smooth host halo density. The total mass of the clump is then 
determined by the mass inside R c . 
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Fig. 1. Anatomy of a WIMP clump. Radial distances are drawn on a 
logarithmic scale, with the center chosen to be 10 13 cm. The NFW 
clump (left) shows r c inside R c i. The SIS clump (right) shows, from 
inner to outer, i? m i n , R1/2, and R c a- Both clumps are 10 8 Mo, with 
R cl ~ 1 kpc. 



In order to study the range of possible clump radial pro- 
files, we also consider the case of clumps described by singu- 
lar isothermal spheres (SIS). The density profile in this case 
is given by: 



p(r) = m x n [ — 

. ' r. 



(2) 



By contrast with NFW, the SIS profile provides a steeper 
slope, leading to a sharp peak density at the center of the 
clump. The SIS model describes a collection of self-gravitating 
collisionless particles. If dark matter clumps lack baryons, 
then this profile is likely to be more appropriate. The pa- 
rameters r c and tiq are fixed by the total clump mass and the 
halo density at the edge. Figure 1 sketches both NFW and 
SIS clumps (i? m in and R1/2 are described in section 4). 

In addition to the specific halo and clump profiles, we need 
a model for the distribution of clumps in the halo. Following 
Blasi and Sheth (2000), clumps of mass m and position r are 
distributed according to 



n c i(r, m) = n c i, 
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where n c i,o is a normalization constant and rf is the core of 
the clump distribution. A value of a ~ 1.9 fits well the sim- 
ulations in Ghigna et al. (1998), in which a halo with Mh ~ 
2 x 10 12 M Q contains about 500 clumps with mass larger 
than M cl min ~ 10 8 M Q . The presence of such dark sub- 
structure gains some support from recent analyses of strong 
gravitational lens systems (see, for example, Chiba (2001)), 
where dwarf galaxy satellites and globular clusters are insuf- 
ficient to affect the lensing, but dark subhalo structures do 
alter the flux ratios in multiply imaged systems. 
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The halo representations generated here extend the clump 
distribution to a lower mass cutoff of 10 7 Mq. In princi- 
ple, lower mass clumps will be present in simulations with 
larger dynamical range. For the case of gamma-ray emis- 
sion, Calcaneo-Roldan and Moore (2001) find that little ad- 
ditional flux arises if the minimum mass cutoff is lowered. A 
more detailed study of lower mass clumps will be described 
in Aloisio et al. (2002b). 

2.2 Generating Synchrotron Radiation 

Neutralinos in DM clumps annihilate mostly into quark-anti- 
quark pairs which hadronize mostly into pions. The neutral 
pions give rise to the gamma-ray emission via ir° — > 77; 
while charged pions decay into charged leptons: 



(4) 



and muons subsequently decay into electrons and positrons: 



e z/„z/ e ,fj, 



(5) 



Neutrinos from this cascade process may eventually be de- 
tectable, while electrons and positrons can be observed via 
synchrotron radiation in the presence of magnetic fields. 

The number of electrons and positrons produced at each 
energy resulting from a single \X annihilation, j^, can be 
calculated via: 
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This form for the pion spectrum is the result of the fragmen- 
tation and hadronization of partons, as found in Hill (1983). 
More sophisticated fits can be used but for the present work 
eq. (7) is sufficient. It is worth noting the strongly non- 
thermal form of the fragmentation spectrum, which is in- 
strumental in distinguishing DM annihilation from the back- 
ground provided by the CMB radiation. In the integral (6), 
the pion decay generates the following muon number per unit 
energy: 
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and the muon decay gives 
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At low energies, eq. (6) behaves as 

dN e 
dE P . 



e: 



(8) 



(9) 



(10) 



and goes to zero as E e reaches the maximum energy which 



is of the order of m x . 



The total injection rate of electrons and positrons by neut- 
ralino annihilations is then given by 



q e (E e ) = n 2 x {av) xx 



dNe 

dE P . 



(ID 



where (av) xx is the annihilation rate. Depending on the spe- 
cific neutralino model, the annihilation rate varies around a 
fiducial value of (av) xx ~ 3 x 10~ 29 cm 3 /s, with a wide 
spread. The electron-positron population at any point in space 
can be written as 



dn e 
~d~Ee 



1eT, 



(12) 



where r is the average lifetime of the electron or positron. 
For the cases considered below, the electron or positron gy- 
roradius is usually much smaller than the size of the emitting 
region. It is a good assumption to calculate the emission of 
electrons and positrons as coming from their original posi- 
tion, t is the timescale over which electrons and positrons 
radiate away most of their energy. 

Depending on the density profile of the clump, the electron 
or positron energy, E e , and the position within a clump, t is 
determined by the shortest among the following timescales: 
the energy loss timescale for (1) synchrotron radiation, (2) 
inverse Compton scattering (ICS) off the cosmic microwave 
background (CMB), (3) ICS off the local synchrotron pho- 
tons, and (4) the possibility of electron-positron annihilation. 
In the NFW case discussed below, the synchrotron emission 
is the main energy loss process. The SIS case has different 
regimes with the central regions being dominated by ICS off 
the generated synchrotron photons. We discuss each case in 
more detail below. 

The total synchrotron power (Rybicki and Lightman, 1979), 



dE s . 



dt 



= 4 x 10" 



l BlE 2 



erg/s, 



(13) 



where is the magnetic field strength in /iGauss and E e is 
the electron energy in GeV. In our Galaxy, estimates of 
vary between 3 and 6. If we choose E e ~ m x ^100 GeV, 
the typical synchrotron power becomes 4 x 10~ 17 erg/s. The 
peak frequency, 



= 3.7 x 10 b B u Et Hz . 



(14) 



in the Galactic /iGauss fileds, indicate that the synchrotron 
emission should be observable below about 50-100 GHz. 

In the numerical calculations of individual clump signa- 
tures discussed below, we model the Galactic magnetic field 
B(r, z) using Stanev (1997), where here r and z are Galaxy- 
centered cylindrical coordinates. We also assume for sim- 
plicity that the synchrotron radiation is mainly produced at 
the peak frequency defined above. 

If the only significant loss process is the emission of synch- 
rotron radiation, then 



Ee 



dE syn /dt ' 



(15) 
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The synchrotron emissivity from annihilating neutralinos is 
then 

. ^ dn e dE e dE syn erg 
dE e dv dt cm 3 s Hz 
In some regions of DM clumps, other processes affect the 
electron-positron number density per energy more strongly 
than the synchrotrom emission. In the SIS case, the most 
significant process is inverse Compton scattering (ICS) off 
the generated synchrotron photons. Other effects come from 
the e+ — e~ pair annihilation and synchrotron self absorption. 
We discuss these as they become relevant in section 4. 

2.3 Comparing with the CMB 

As discussed above, the range of frequencies where synch- 
rotron emission from DM clumps may be observed peaks 
around <~ 40 GHz (m x /100GeV) 2 . This frequency range is 
the focus of intensive studies of the structure of cosmic mi- 
crowave background (CMB) radiation. The CMB signals are 
well understood and well measured and give the opportunity 
to search for the DM clump emission. We focus our atten- 
tion on the band between 10 — 400 GHz, which should be 
relatively free of contaminating CMB foregrounds (Tegmark 
et al., 2000), and which is where the most sensitive CMB 
experiments are planned to operate. 

Only highly peaked SIS clumps may reach the levels of the 
CMB emission. In general, the annihilation signal is compa- 
rable in flux to the CMB anisotropies, which are at the level 
of <~ 10~ 5 of the CMB. In what follows, we assume that the 
CMB anisotropies have a blackbody spectrum. 

The best spectral data on the CMB come from the COBE 
FIR AS instrument (Fixsen et al., 1996). The reported uncer- 
tainties are of order of 20 kJy sr~\ above 68 GHz and only 
above 7° angular scales. The COBE DMR instrument (Smoot 
et al., 1992) used 31, 53, and 90 GHz, and found consistent 
anisotropies (of higher order than dipole) at the 1.1 x 10~ 5 
level, but again on 7° scales. A combined FIRAS and DMR 
analysis (Fixsen et al., 1997) sets the anisotropy spectrum at 
about 5 x 10~ 5 of the cosmic monopole spectrum for these 
large scales. 

Our interest is in pixel sizes typical of upcoming CMB 
probes. The MAP mission (Hinshaw, 2000) is typically 18 ar- 
cmin and the Planck Surveyor (de Zotti et al., 1999) is 5 ar- 
cmin, so we choose a nominal pixel size of 10 x 10 arcmin 
(the smaller the pixel, the brighter the clump, because one 
can aim at the central cusp where the neutralino annihilations 
are most rapid). The BOOMERANG (de Bernardis et al., 
2000) and MAXIMA- 1 (Hanany et al., 2000) experiments 
employed this resolution, and they report 2.6 x 10~ 5 fluctu- 
ations at 1° (BOOMERANG) and 1.7 x 10~ 5 at 14.4 arcmin 
(MAXIMA- 1), with both quoted at 150 GHz. 

On the lower frequency side, some of the best measure- 
ments come from the VLA, probing very small angular scales 
(where primordial plasma fluctuations become damped dur- 
ing the nonzero timescale of recombination). They report 
2 x 10~ 5 at 8.4 GHz (Partridge et al., 1997; Fomalont et al., 
1993). 



The VLA measurements compare meaningfully with clumps 
that appear as point sources. For extended sources, is seems 
best to blend the various experimental results, keeping in 
mind the pixel size of interest. Fortunately, for all scales, 
the data are consistent to within an order of magnitude. We 
choose A/„//„ ~ 2 x 10~ 5 as the limit for detection of de- 
viations distinct from those inherent in the CMB, with the 
intended caveat that actual detections made in this way de- 
serve further scrutiny regarding the precision of the CMB 
blackbody. 

Then to summarize, our criteria for clump detection over 
the CMB are as follows: Between 10 and 400 GHz, the cen- 
ter pixel of the clump (taken to be 10 x 10 arcmin) must 
exceed 2 x 10~ 5 times J^cmb in the same pixel solid angle. 

3 NFW Clumps 

The NFW density profile, given in eq. 1, has a lower density 
central cusp when compared with a SIS clump. This makes 
the NFW case simpler: absorption of synchrotron photons by 
relativistic electrons is unimportant and e pair annihilation 
occurs on a longer time scale than synchrotron radiation, so 
it can be ignored as well. We also neglect the presence of any 
constant density core, as this region is in general too small to 
be detected. 

In order to derive the NFW clump signature, it is sufficient 
to evaluate the integral for the flux observed on Earth, from 
a clump at distance d: 

1 /"^ cl erff 

T - = T^2 4 ^ r ^ dr 2 XT ' < 17 > 

47rcr ,/ crrr s Hz 

using eq. 16. Because r c is relatively large (^75 pc for a 
10 8 M Q clump), NFW clumps will often have observable 
angular sizes. 

The detection of a NFW clump can be made by observ- 
ing both the spectrum of the synchrotron emission as well 
as the detail angular profile which is resolvable for NFW 
clumps. In figure 2, we plot the spectrum of a specific NFW 
clump with mass 10 8 M Q located at Galactocentric coordi- 
nates [-4,0,0] kpc (the Sun is at [-8.5,0,0]), and with the fol- 
lowing base properties: (<jv) xx = 3 x 10~ 29 cm 3 /s and 
m x = 100 GeV. The thick line gives this clump's spectrum 
from its central pixel, and the thin solid line shows the CMB 
anisotropies (taken as 2 x IO^J^cmb) in the same pixel 
solid angle. The other lines show the effect of modifying the 
mass of the clump by a factor of 10 (dashed line), increas- 
ing the cross section by a factor of 100 (dotted line), and 
changing the neutralino mass by a factor of 10 (dot-dashed 
line) or a factor of 100 (dot-dot-dashed line). One can see 
that raising m x results in a higher cutoff frequency, so that 
if a bright enough annihilating neutralino clump were ob- 
served, its spectral cutoff would immediately reveal m x . On 
the other hand, changes to M c \ and (<Jv) xx are degenerate, 
so some other means is necessary to distinguish between the 
two, such as gravitational lensing, or dynamical effects on 
nearby stars. 
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Fig. 2. Spectrum of synchrotron emission from an example NFW 
clump, compared with the CMB anisotropies (thin line) in the same 
solid angle as the clump. This clump is chosen to lie at Galacto- 
centric coordinates [-4,0,0] kpc, 4.5 kpc away from us. The heavy 
solid line is the base case, where M c \ = 10 8 Mq, (cv) X x = 
3 x 10 -29 cm 3 /s, and m x = 100 GeV. The dashed line increases 
M c i by 10 times; the dotted line increases (<jv) xx by 100 times; the 
dot-dashed line increases m x by 10 times; the dot-dot-dashed line 
increases m x by 100 times. 



We have created several simulations of a clumpy halo. In 
each, the lower cutoff on the clump mass, M c i m i n , is a spec- 
ified parameter. Table 1 presents the major results for each 
of six halo realizations with M c i lTO i n = 10 7 Mq, for which 
3972 clumps were generated, totaling ~ few 10 11 Mq, These 
halo simulations all employ base case attributes ((crv) xx = 
3 x KT 29 cm 3 /s and m x = 100 GeV). 

We find that between 250 and 300 clumps will be observ- 
able in this halo, above the CMB anisotropies between 10 
and 400 GHz as described in section 2.3. None of the NFW 
clumps is bright enough to outshine the CMB radiation it- 
self; only the CMB fluctuations are exceeded. The NFW 
clumps in these realizations occupy slightly less than 0.1% 
of the full sky's solid angle (assuming a minimum pixel size 
of 10 x 10 arcmin), and typically several clumps are within 
3 kpc of the Earth. (Much of this information in table 1 is 
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Fig. 3. A simulated sky of observable NFW dark matter clumps. 
Each cross indicates a clump, observable by the criteria laid out 
in section 2.3. The size of the cross indicates the angular size of 
the clump such that half its light originates inside the area shown. 
Angular sizes are exaggerated by a factor of 10. 



intended for comparison with the corresponding SIS case, 
given in table 2 and discussed in section 4.3.) 

Of the few hundred observable clumps, between ~ 50 and 
75 will reside above 30° or below —30° Galactic latitude, 
where Galactic microwave contamination is not problematic. 
The locations of all the clumps in the sky are depicted for 
halo realization 1 in figure 3, with crosses indicating the an- 
gular size of the clumps such that half of their total emis- 
sion comes from inside that area. This half-light angle cor- 
responds to a half-light radius, which we will call R\/2, that 
equals 0.26r c for an NFW clump. For clarity, the angular 
sizes in the figure are exaggerated by a factor of 10. 

A histogram of the number of clumps vs. angular size is 
given in figure 4. The figure shows results for all 6 halo real- 
izations discussed so far, plus several other cases correspond- 
ing to modified parameters M c \^ n , (av) xx , and m x ; we 
modify these in order to cover some of the possible halos and 
some of the supersymmetric parameter space available to the 
neutralino. In each case, ~ 100 clumps are resolvable, un- 
less m x is diminished by ~ 100 times which is too small for 
any reasonable neutralino model. Interestingly, if the clump 
mass spectrum extends below 10 7 Mq (in which case there 
are 23,830 clumps still totaling — few x 10 11 Mq), or if 
(o~v) xx is significantly stronger than assumed, many more 
clumps can be found. A angular correlation function for this 
case is studied in Aloisio et al. (2002b). 



Table 1. Some results for NFW clumps in six halo realizations, 
each set with M c \ >m i D — 10 7 Mq, yielding 3972 clumps. The 
following information on observable clumps is given: number of 
them, number which outshine the CMB (rather than just the CMB 
anisotropies), number above 30° or below —30° Galactic latitude, 
number within 3 kpc of Earth, and the percentage of the sky the 
clumps occupy. 



4 SIS Clumps 

The SIS density profile was given in eq. 2. The steep cusp 
diverges at the center and makes it necessary to define i? m in, 
the radius inside which neutralino annihilations are so rapid 
that n x (r < R m i n ) = n x (R m i n ) remains constant. i? m i n 
is found by setting the cusp-forming timescale equal to the 
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Fig. 4. Histogram showing the half-light angular size of NFW 
clumps in several halo realizations, as used in figure 3. Solid lines 
are for halo realizations 1 — 6, where M C L m i n = 10 M©. The 
dashed line is for M c i jm j n = 10 Mq; the dotted line is for {<rv) xx 
increased 100 times; the dot-dashed and dot-dot-dashed lines are 
for m x increased by 10 and 100 times, respectively. The thin dotted 
vertical line shows the angular scale of a 10 x 10 arcmin pixel for 
reference. 



XX annihilation timescale (Berezinsky, Gurevich, and Zybin, 
1992), so that 



-Rmin = -Rcl {(TV) 



A X 



^halo 



(18) 



where rihaio is the halo DM density at the location of the 
clump. 

In the absence of other interactions, one would integrate j v 
from eq. 16 over the volume of the clump, and find that the 
dominant part of the flux comes from the vicinity of i? m i n - 
However, due to the high density at the center of an SIS 
clump, there are other important factors to include. We dis- 
cuss each factor in turn: 

4. 1 Relevant Interactions 

(1) Inverse Compton scattering of electrons and positrons 
against the synchrotron photons created by the electrons and 
positrons is the most important correction to the electron- 
positron distribution in space and energy. To see this, we 
write the ICS power as (Rybicki and Lightman, 1979) 



(19) 



1, and 



where <tt = 6.65 x 10" 

(which is appropriate in the limit hu <C m e c 2 ), (3 
t/ph is the energy density in the photon field. 

The photon population with the greatest energy density 
controls the rate of ICS. For the CMB, 
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where a = 7.56 xlO" 15 erg cm~ 3 K~ 4 and T cmb = 2.728 K. 
At the inner cores of the SIS clumps, we find that the synch- 
rotron photon population dominates the photon energy den- 
sity such that J7 p h — U syn > U cm b- The synchrotron energy 
density at distance r from the clump center can be calcu- 
lated as the integral of the emissivity along all lines of sight, 
s(r), from r to all other points in the clump. The energy 
and space dependence in the equilibrium density of elec- 
trons and positrons can be factorized as dn e (E e , r)/dE e oc 
<7o(r) following from a similar factorization in the in- 

jection function q(E e ,r). This leads to 



C/ syn (r) cx B 2 m 3 x I(r) 



where 



1 f 1 f s ^ 

I(r) = — / dcosO I a {r{s))ds 

2 J-i Jo 



(21) 



(22) 



and s max = rcosO + {{rcosO) 2 + R 2 - r 2 ) 1/2 . After some 
algebraic steps we find that 



<7 (r) oc r /I(r), 



(23) 



where j3 is the slope of the radial dependence in the injection 
function ((3 = 4 for isothermal clumps outside the minimum 
radius). The detailed equation is then solved iteratively in 
order to calculate the equilibrium electron positron function, 
which gives 



dn e (E ei r) 
dE K 



oc E, 



-5/2 r -5/2 



(24) 



for an isothermal clump. Note that in the most general case 
that we considered numerically, the equation to solve is an in- 
tegral equation since the function J(r) contains the unknown 
equilibrium distribution. 

This holds wherever electron losses are fastest by ICS against 
the synchrotron photon field. In the case of losses being dom- 
inated by synchrotron radiation and ICS against the CMB in- 
stead, the loss mechanism is independent of r, and dn e /E e oc 



E, 



-5/2 



Since the two forms of electron losses have dif- 



ferent power laws in r (and the same power law in E e ), there 
is a transition radius i?i cs ,syn (which is independent of E e ) 
between the two behaviors. 

(2) Pair annihilation between e~ and e+ particles would 
have an important role in determining the flux from a DM 
source in the absence of the ICS described above. But in the 
presence of ICS losses, pair annihilation makes little differ- 
ence to the resulting radio flux. Pair annihilation is still im- 
portant, because it reduces the population of the lower energy 
pairs, which would otherwise upscatter photons frequently 
enough to affect the synchrotron spectrum. 

We can estimate the e ± lifetime for the case where elec- 
trons annihilate before they lose a significant fraction of their 
original energy via synchrotron or ICS as 



T p + 



{n e 



1 



(25) 
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with v e ~ c, and the angle brackets indicating an average 
over E e ~. (This is equally valid if we interchange all the 
plus and minus signs, but we will keep track of signs in order 
to distinguish between particles and antiparticles.) 
With eq. (12), we have 

dn e + = g e + (26 ^ 
dE e + (n e -a e ±v e -)' 

The term in angle brackets is computed from the equations 
in Svennson (1982) (also consistent with Coppi and Bland- 
ford (1990)). (n e -a e ±v e -) can be expressed as an energy 
integral which depends on the e ± distribution that we are 
trying to determine. Again, we find the solution to this non- 
linear equation via numerical iterations. 

Neglecting a slow varying logarithm in the annihilation 
cross section, one finds that n e - is roughly constant and 
(<r e ±v e -) oc E~ + . That is, a given positron of any energy 
sees essentially the same distribution of target electrons, so 
to a first approximation, it is sufficient to replace the elec- 
tron distribution with a single population at energy E e - ~ 
m e c 2 . By that rationale, r e + ~ E e +. At low energies where 
dN e /dE e ~ Ee 3/2 , this means that dn e /dE e ~ E^ 1 / 2 and 
n e ~ E 1 / 2 (where the subscript e refers to either an electron 
or a positron). 

The distribution in space, rather than energy, is easier to 
determine. In order to satisfy eq. (26) with the injection 
q e oc r~ 4 set by the isothermal DM profile, we must have 
dn e /dE e oc r~ 2 . Then for a region whose e ± losses are 
fastest by pair annihilation, we have 

5"^' <27 > 

We define i? arm .ics as the transition radius between the 
r~ 2 pair annihilation distribution and the r~ 5/>2 ICS distri- 
bution. If i?ann,ics > -R m i n , pair annihilation dominates the 
central region over ICS. In our halo and for most choices 
of m x , i?ann,ics < ^?min for electrons radiating at frequen- 
cies of interest. The distribution flattens at the center before 
the density can get high enough to create a pair annihilation 
dominated region. At lower frequencies, pair annihilation be- 
comes important and limits the number of low energy pairs 
which become unable to effect the synchrotron spectrum. 

(3) Synchrotron self absorption (SSA) occurs when the 
electrons in the magnetic field can absorb the energy of the 
synchrotron photons. The per-unit-length absorption coeffi- 
cient for SSA can be written as: 



dE, 



Op 



8m/ dt 



d_ 

dv 



1 dn e 
v dv 



(28) 



where the photon frequency is related to the electron energy 
via v max . SSA absorbers of a particular photon are electrons 
of similar energy to the original synchrotron emitter which 
made the photon. As such, SSA is only effective below some 
^ cr it, where the corresponding electrons are more numerous. 

The source function for an absorbing source is S v — j v /a v , 
and the optical depth is 




10° 10" 10 

radial coordinate (cm) 

Fig. 5. Radial flux profile for an example SIS clump, integrated 
along the line of sight through the clump. This plot gives the relative 
intensities found by pointing toward different parts of the clump, 
given sufficient angular resolution. The solid line is for 10 GHz and 
the dashed line is for 100 GHz. 



where z is a line-of-sight coordinate. The flux density from 
such a source at distance d from the Earth is obtained by 
summing over each line of sight through the DM clump: 



47TC2 2 Jo 



2nbS v {l 



T ")db 



erg 



cm 2 s Hz 



(30) 



The SSA cutoff frequency v crlt depends in part on the ab- 
sorber density, and without ICS or pair annihilation ef- 
fects, SSA would be the most important correction to the 
calculated flux. SSA has been included in the calculations 
performed here, but since the central absorber population is 
diminished by ICS, it is of secondary importance. 

4.2 Results 

The synchrotron spectrum of DM clumps is predominantly 
altered by ICS off synchrotron photons. Figure 5 shows an 
example of a SIS clump profile that shows the flattened cen- 
tral region. A characteristic size for the clump, R\/2, can 
be defined by the radius inside which half of the synchrotron 
radiation originates. 

Figure 6 shows the spectrum of synchrotron emission from 
a SIS clump in the same position as the NFW clump of fig- 
ure 2. The thick solid line gives the base case, M c \ — 10 8 M©, 



dz , 



(av) xx = 3 x 10~ 29 cm d /s, and m x = 100 GcV. The 
dashed line increases M c \ by 10; the dotted line increases 
(<jv} xx by 100 times; the dot-dashed line increases m x by 
10 times; the dot-dot-dashed line increases m x by 100 times. 
The thin solid line gives the intensity of the CMB within the 
10 x 10 arcmin pixel, for comparison. 

In order to compare the SIS with the NFW clumps, we 
used the same halo simulations as in the NFW case: realiza- 
tions 1-6 with 3972 clumps (M c i jmin = 10 7 M ), and an- 
other realization with 23,830 clumps (M c i im ; n = 10 6 M©), 
(29) a n w ith ~ few x 10 11 M© total. For a typical clump mass of 
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10 10 
frequency (GHz) 



Fig. 6. Spectrum of synchrotron emission from an example SIS 
clump chosen to lie at Galactocentric coordinates [-4,0,0] kpc. The 
thick solid line gives the base case, M c \ = 10 8 Mq, (o~v) xx = 



3 x 10" 



/s, and 



= 100 GeV. The dashed line increases 



Md by 10 times; the dotted line increases ■ 



by 100 times; the 



dot-dashed line increases m x by 10 times; the dot-dot-dashed line 
increases m x by 100 times. The thin solid line gives the intensity 
of the CMB within the 10 x 10 arcmin pixel, for comparison. 



10 8 Mq, a SIS clump can outshine a NFW clump by about 
10 5 times in the microwave. 

The visible SIS clumps for halo realization 1 are shown in 
figure 7, here with angular sizes based on R\/2 exaggerated 
by a factor of 5,000 (none of these clumps is resolvable with- 
out interferometry). In table 2, we list some results of the SIS 
case. First, we find that SIS clumps (in an NFW halo) tend to 
produce 600 — 700 observable sources, which is a few times 
more than the NFW clump case (compare tables 1 and 2). 
The combined sources occupy only slightly less solid angle 
than with NFW, mostly due to the small angular size of SIS 
clumps. But unlike the NFW case, ^20% of the SIS clumps 
not only outshine the CMB anisotropics, but outshine the 
CMB altogether. A similar percentage of observable clumps 
is located outside the direction of the Galactic disk. Finally, 
the table gives several observable clumps within 3 kpc of the 
Earth, which can be resolved (as in figure 5) by long baseline 
interferometers. 

Table 3 gives a comparison of the NFW and SIS clump 
profiles, by reporting the number of observable clumps in 
each case, while also varying m x , (av) xx , or M c i, m j n . The 
results show that changing these parameters has different out- 
comes for NFW and SIS. For example, increasing m x by a 
factor of 100 lowers the spectrum for most relevant v bands 
(see figure 2); therefore, NFW clumps get dimmer and fewer 
are seen. On the other hand, the SIS clumps are bright enough 
to avoid that worry, and the increase in m x also increases the 
maximum electron energy. This allows strong sources even 
in weak magnetic fields via eqs. 13 and 14, raising the num- 
ber of observable clumps. 

As another example, changing {vv) xx has little effect on 
SIS clumps, because they have an overdensity of electrons 



+■ 



■*, 



-60 -30 3 60 
Galactic azimuth (deg) 



120 150 180 



Fig. 7. A simulated sky of observable SIS dark matter clumps. Each 
mark indicates a clump, observable by the criteria laid out in section 
2.3. The size of the cross indicates the angular size of the clump 
such that half its light originates inside the area shown. Angular 
sizes are exaggerated by a factor of 5000. 



and positrons emitting anyway. But for NFW, where the 
clumps are typically near the CMB anisotropy level, chang- 
ing the XX cross section affects the fluxes significantly. Both 
NFW and SIS cases see similar growth in the number of ob- 
servable clumps when the clump mass spectrum is extended 
downward. 



5 Conclusion 

In this paper, we studied the synchrotron signature of neut- 
ralino annihilation in the presence of the Galactic magnetic 
field. If hierarchical CDM clustering extends down to mass 
scales smaller than dwarf galaxies, these DM clumps are de- 
tectable in microwaves over the fluctuations in the CMB. 

If the clumps have nearly NFW profiles like the Galaxy, 
then the clumps have finite angular sizes. NFW clumps have 



SIS 


1 


2 


3 


4 


5 


6 


observable 


636 


632 


676 


656 


682 


616 


outshine CMB 


54 


34 


38 


37 


37 


40 


> 30° latitude 


130 


155 


144 


129 


140 


135 


within 3 kpc 


9 


8 


2 


5 


5 


7 


resolvable 


212 


184 


209 


190 


210 


201 


resolve > 30° 


48 


53 


51 


49 


54 


59 


% of sky 


0.04 


0.04 


0.05 


0.04 


0.05 


0.04 



Table 2. Some results for SIS clumps in six halo realizations, 
each set with M c i jm i n = 10 7 Mq, yielding 3972 clumps. The 
following information on observable clumps is given: number of 
them, number which outshine the CMB (rather than just the CMB 
anisotropies), number above 30° (or below —30°) Galactic latitude, 
number within 3 kpc of Earth, number resolvable by 10 km interfer- 
ometry such as the VLA, number resolvable by the VLA above 30° 
Galactic latitude, and the percentage of the sky the clumps occupy. 
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case 


NFW 


SIS 


base case 


282 


636 


(av) xx X 100 


528 


671 


m x x 10 


284 


1126 


m x x 100 


66 


917 


M c i !min x 0.1 


1078 


3503 



Table 3. Comparison of NFW and SIS clumps. The base case refers 
to m x = 100 GeV, (av) xx ~ 3 x 10~ 29 cm 3 /s, and Af c i, min = 
10 7 Mq . The other cases modify these parameters as stated. The 
last row of the table uses a different halo realization than the other 
rows. 

synchrotron spectrum from non-thermal electrons with neg- 
ligible ICS and self-absortion. NFW clumps will be large 
enough and numerous enough that several should be resolv- 
able in angular profile by CMB anisotropy experiments. The 
light profile will follow r _1 inside a core radius. The statis- 
tics of the clumps distribution should be consistency with fig- 
ure 4 and table 1 . 

SIS clumps are brighter and smaller in angular size. The 
clumps have a distinctive spectrum due to a combination of 
ICS and pair annihilation features, in a particular radial struc- 
ture. The spectrum can be checked at lower frequencies also, 
because these clumps are very bright and may outshine other 
radio backgrounds. A few clumps should be close enough to 
display a distinctive radial profile by interferometric means. 

In either case, one should make use of CMB anisotropy 
studies to search for neutralino annihilation in DM clumps. 
Many clumps should have signal above 2 x 10~ 5 times the 
2.728 K Planck function. The sources are usually small and 
occupy a small fraction of the sky, so that a full sky high res- 
olution microwave survey maybe necessary. The upcoming 
MAP and Planck missions will provide that. Even in smaller 
data sets such as the BOOMERANG, MAXIMA, and DASI, 
it is possible that clumps can be detected. 

Upon finding synchrotron source clumps or hot pixels, one 
final check can prove their DM nature. The gamma ray and 
neutrino emission from them must fit the expected levels, 
and must be independent of position in the Galaxy, unlike 
the synchrotron which follows the Galaxy's magnetic field. 
This combination of signatures makes the sources unique, 
and aids in extracting specific information about them: M c \, 
m x , and (av) xx . 

One interesting point to note is that for reasonable choices 
of the appropriate parameters, we find that several hundred or 
more clumps are typically observable in either profile. If new 
CMB experiments fail to detect them, then we have learned 
that either neutralinos are not the dark matter, or CDM halo 
clumps do not exist. However, if they are found, they offer 
a valuable opportunity to learn more about CDM clustering, 
neutralino properties like mass and cross section, and possi- 
bly the Galactic magnetic field structure. 
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